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METHODOLOGY  FOR  ESTIMATING  SMOKE/QBSCURANT 
MUNITION  EXPENDITURE  REQUIREMENTS 


I.  INTRODUCTION. 

For  the  last  few  years,  the  US  Amiy  has  expended  a large  effort  to  improve  its  ability 
to  create  and  utilize  smoke/obscurant  screens  tn  the  battlefield  against  enemy  electro-optical  target 
acquisition  devices.  A key  aspect  in  the  research  and  development  process  is  the  analysis  of 
potential  effectiveness  of  developmental  systems.  Since  most  of  the  technology  of  electro-optics  has 
been  developed  since  World  War  II  and  since  little  interest  was  paid  to  smokes/obscurants  until  the 
1973  Arab-Israeli  conflict,  the  methodology  to  support  systems  analysis  studies  has  been  developed 
only  recently.  This  methodology  has  been  concentiated  in  the  many  models  which  attempt  to 
simulate  smokes/obscurants  in  the  battlefield.  Techniques  have  been  developed  to  relate  simulated 
performance  to  measures  of  effectiveness,  but  little  documentation  of  these  techniques  has 
appeared.  The  purpose  of  this  report  is  to  disseminate  some  analytical  tools  developed  at  Chemical 
Systems  Laboratory  to  support  systems  analysis  studies  of  developmental  smoke  munitions. 

Smoke  munitions  effectiveness  estimates  are  usually  based  upon  munition  expenditure 
requirements  to  accomplish  a screening  mission.  The  mission  is  normally  described  in  terms  of  the 
target  to  be  screened  from  observers  located  along  a front;  i.e.,  the  specific  location  of  the 
observers)  is  assumed  to  be  unknown,  therefore  a front  line  must  be  screened.  The  target’s 
background,  the  natural  and  artificial  sources  of  radiant  illumination  may  be  specified,  along  with 
the  target  acquisition  device,  e.g.,  unaided  eye,  binoculars,  electro-optical  (EO)  imager.  A variety  of 
meteorological  conditions,  which  strongly  affect  munition  performance,  may  be  considered.  In  an 
idealized  study,  only  crosswind  and  headwind/tailwind  conditions  are  considered.  Quartering/ flank 
wind  conditions  lead  to  expenditure  requirements  between  the  two  extremes. 

For  a specified  screening  mission,  there  is  probably  a minimal  set  of  optimum  aiming 
points  at  which  multiple  rounds  should  be  fired  to  accomplish  the  mission.  However,  the 
complexity  of  simply  stating  the  associated  constrained  optimization  problem  is  beyond  practical 
consideration.  One  might  consider  a trial-and-error  approach  by  simulating  the  firing  of  multiple 
rounds  at  a set  of  arbitrary  aiming  points.  A more  practical  approach  is  to  estimate  the  performance 
of  a single  munition  (or  submunition)  under  the  conditions  of  the  mission  and  use  some  systematic 
method  to  estimate  the  number  of  munitions  required  to  screen  the  front.  The  development  of  such 
methods  is  the  subject  of  this  report. 

II.  BASIC  CONCEPTS. 

A.  Instantaneous  Sources. 

When  a smoke  agent  is  released  over  a very  short  time,  say  less  than  1 second,  it  is 
considered  an  instantaneous  volume  source  with  some  initial  cloud  dimensions.  If  a horizontal  plane 
containing  all  lines  of  sight  at  a given  height  passes  through  the  cloud,  the  intersection  of  the  cloud 
boundary  with  the  plane  is  a closed  curve  bounding  an  area.  For  a Gaussian  cloud,  which  is  assumed 
in  this  methodology,  the  effective  cloud  boundary  is  defined  by  the  set  of  points  at  which  some 
property  of  the  cloud  has  a specified  minimum  value.  Depending  upon  the  obscuring  mechanism 
involved,  the  critical  property  may  be  concentration,  the  integrated  concentration  (CL),  the 
reduced  contrast  resulting  from  the  obscuring  power  of  the  cloud,  or  the  transmittance  through  the 


cloud.  The  maximum  width  of  the  bounded  area,  measured  normal  to  the  line  of  sight,  is  defined  as 
the  effective  screen  width  of  the  cloud  at  the  height  of  the  line  of  sight  at  the  instant  in  time. 

As  time  increases,  the  cloud  expands  due  to  turbulent  diffusion  and  is  transported 
downwind  at  the  prevailing  windspeed.  The  effect  of  these  phenomena  is  to  cause  the  effective 
screen  to  grow  while  moving  laterally.  The  estimation  of  munition  expenditure  requirements  thus 
consists  of  two  basic  tasks:  ( 1 ) prediction  of  effective  screen  width  as  a function  of  time  and 
(2)  selection  of  the  number  and  locations  of  a set  of  aiming  points  which  will  be  expected  to  cause  a 
specified  front  to  be  screened  as  of  some  time  after  the  initial  round  or  rounds  function. 

B.  Continuous  Sources. 

The  most  effective  smoke  screening  systems  today  make  use  of  continuous  sources, 
which  emit  smoke  over  a long  time  interval.  Theoretically,  smoke  generators  can  operate 
indefinitely,  but  actually  all  real  systems  emit  for  a finite  time,  usually  referred  to  as  the  burning 
time  (tg).  Many  new  munitions  deliver  a set  of  submunitions,  ranging  in  number  from  6 to  92. 

A single  continuous  source  produces  a smoke  cloud  shaped  like  a long  cone  lying  on  its 
side  with  the  vertex  at  the  source.  The  base  end  of  the  cone  may  at  some  time  slow  down  in  its 
growth  rate  and  tend  to  form  a cylindrical  plume  parallel  to  the  wind  direction.  If  a crosswind  line 
of  sight  is  to  be  screened,  i.e.,  a line  of  observers  parallel  to  the  wind,  the  effective  screen  length  is 
defined  as  the  distance  between  the  farthest  downwind  and  the  least  downwind  lines  of  sight  which 
are  effectively  screened.  For  a headwind/tailwind  line  of  sight,  the  effective  screen  width  is  defined 
in  the  same  way  as  for  an  instantaneous  source.  The  expected  variation  with  time  is  different, 
however,  and  allows  for  a simpler  technique  to  estimate  munition  expenditure  requirements. 

For  munitions  with  multiple  submunitions,  the  effective  crosswind  screen  length  is 
defined  the  same  as  for  a single  source,  but  there  is  no  practical  way  to  combine  the  effects  of 
multiple  clouds  from  randomly  located  sources  by  other  than  simulation  techniques.  The  normal 
approach  to  the  simulation  would  be  a Monte  Carlo  technique  varying  the  orientation  of  the 
elliptical  pattern  of  submunition  locations  (impact  points)  to  the  wind  direction.  Hundreds, 
possibly  thousands,  of  effects  estimates  would  be  averaged  to  obtain  an  expected  effective  screen 
length.  Some  smoke  simulation  models,  particularly  the  one  used  by  the  author,  are 
computationally  too  expensive  for  Monte  Carlo  methods.  The  procedure  we  use  attempts  to 
estimate  the  expected  set  of  impact  points  which  in  a single  simulation  produce  the  expected 
effective  screen  length.  The  procedure  and  rationale  are  described  below. 

Assume  field  test  data  are  available  for  the  firing  of  a single  round  which  releases 
multiple  submunitions,  each  of  which  is  a continuous  source.  The  impact  points  are  uniformly  and 
randomly  distributed  in  an  elliptical  area.  The  ellipse  has  semi-axes  of  lengths  a and  b.  The  area  of 
the  ellipse  is  7rab.  The  rotation  of  the  ellipse  about  its  center  would  cause  the  ends  of  the  axes  to 
circumscribe  concentric  circles  of  radii  a and  b.  These  circles  describe  the  minimum  and  maximum 
dimensions  of  impact  areas  which  a crosswind  line  of  sight  would  “see,”  because  it  will  intersect 
every  downwind  plume  regardless  of  the  crosswind  dimension.  We  take  as  an  estimate  of  the 
expected  impact  area  a circle  whose  area  is  that  of  the  ellipse;  i.e.,  7rr2  = 7rab.  Solving  for  r gives  the 
radius  of  the  circle  as  the  geometric  mean  of  the  semi-axes;  i.e.,  r = ^/ab.  Using  this  circle,  one  can 
develop  a technique  for  generating  a set  of  impact  points,  relative  to  the  center  of  the  circle,  which 
are  uniformly,  but  not  randomly,  distributed.  Once  a characteristic  impact  pattern  is  generated  for  a 
multiple  submunition  round,  the  determination  of  effective  screen  length  is  the  same  for  all  systems 
of  continuous  sources. 


111.  METHODOLOGY. 


A.  Instantaneous  Sources. 

I . Effective  Screen  Width. 

The  task  of  predicting  effective  screen  widths  with  time  may  be  accomplished  in  many 
ways  because  there  are  many  different  models  available  to  simulate  smoke  munition  performance. 
The  reader  may  have  his  own  technique,  but  a procedure  using  ACT  (ASL/CSL/TRASANA),  a 
derivative  of  the  original  SOM  I developed  by  Johnson,*  will  be  described  here. 

An  instantaneous  smoke  source  released  at  the  origin  of  a reference  system  in  which  x 
is  the  downwind  direction,  y is  the  crosswind,  and  z is  the  height  above  ground  is  centered  at  time  t 
at  the  point  [ (ut,  0,  z(t)J.  Its  dimensions  are  expressed  in  terms  of  the  standard  deviations  of  the 
spatial  distribution  of  the  mass:  ox(t),  Oy(t),  and  oz(t).  Since  we  are  concerned  with  horizontal 
screen  width,  az( t)  is  not  involved  in  this  discussion.  The  boundary  of  the  cloud  is  an  ellipse  with 
the  major  and  minor  axes  being  proportional  to  ox  and  oy.  For  the  purpose  of  this  method,  we 
prefer  to  assume  ax(t)  = cry(t),  which  is  usually  not  in  serious  conflict  with  reality  but  has  been 
shown  to  be  false.  The  rationale  for  this  assumption  is  the  practical  desire  to  use  one  set  of 
computations  for  all  wmd/line-of-sight  orientations.  With  ax(t)  = ay(t),  the  cloud  has  a circular 
cross  section  which  yields  the  same  effective  screen  width  for  all  lines  of  sight. 

Let  us  denote  by  H(t)  the  half-width  of  the  effective  screen.  If  one  exercises  ACT  with 
a set  of  lines  of  sight  parallel  to  the  x-axis  and  for  a suitable  set  of  t values,  assuming  for  this 
illustration  that  relative  contrast  (C)  is  the  critical  parameter  to  effectiveness,  then  a matrix  of 
contrast  values  can  be  obtained  which  has  the  form: 
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’Johnson,  Morris  C.,  and  Forney,  Paul.  Effectiveness  of  Obscuring  Smokes.  MUCOM  Operations  Research  Group 
Report.  August  1972.  Unpublished. 
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Suppose  an  effective  screen  exists  if  C < CEFF.  Then  for  each  tj,  one  can  find  by  interpolation  the 
value  of  y for  which  C(y,  tj)  = CEFF.  (It  has  been  observed  that  lnC  = a + byisa  good  model  to 
use  for  interpolating  with  relative  contrast.  Other  critical  parameters  may  require  other  models.) 
The  values  of  y(t)  obtained  in  this  manner  are  observed  values  of  H(t)  if  the  cloud  is  initiated  at  the 
origin  and  the  yj  > 0 for  all  j = 1,  2,  . . . . , m.  It  is  possible  then  to  use  least-squares  techniques  to 
obtain  a functional  model  of  H(t),  perhaps  of  the  form: 


H(t)  = 


V 


(1) 


The  values  of  m and  n are  dependent,  of  course,  on  the  number  of  parameters  in  the  H(t)  model  and 
the  desired  resolution. 


2.  Munition  Expenditure  Requirement. 

For  any  specified  screening  mission,  there  may  be  many  ways  to  accomplish  the 
mission.  The  real  world,  however,  imposes  many  constraints  in  terms  of  logistics,  human 
performance  limitations,  and  arbitrary  time  limitations.  To  a certain  degree,  the  methodology 
should  accommodate  such  constraints,  but  practical  considerations  may  require  the  use  of 
simplifying  assumptions.  The  assumptions  underlying  this  methodology  are  as  follows: 

a.  The  mission  is  to  form  an  effective  screen  within  tQ  seconds  of  initial 
smoke  release  along  a front  of  length  FRONT  and  maintain  the  screen  for  ST  seconds. 

b.  All  munitions  are  identical  in  their  performance. 

c.  The  delivery  system  is  capable  of  putting  a round  on  a specified  aiming 
point  (zero  system  error). 

d.  The  fire  unit  has  NTUB  delivery  systems  (artillery  tubes,  mortars. 

launchers). 

e.  Each  delivery  system  may  be  aimed  and  fired  independently  or  in  salvos, 

as  required. 

f.  The  fire  unit  has  a total  load  of  LD  rounds  available  for  the  mission. 

g.  The  minimum  time  between  successive  fires  is  TMU  seconds. 

h.  The  fire  unit  can  maintain  its  maximum  fire  rate  for  TUM  seconds,  after 
which  no  fire  is  possible. 


i.  The  aiming  points  must  be  at  least  100  m from  friendly  troops  to  avoid 
accidental  friendly  casualties  from  impacting  munitions. 
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j.  A fixed  geometry  exists  with  the  lines  of  sight  normal  to  two  parallel 
fronts,  one  for  friendly  targets  and  one  for  enemy  observers. 

k.  The  wind  direction  forms  an  acute  angle  0 with  the  fronts, 

i.e..  0 ^ 0 < 90°. 

These  simplifying  assumptions  lead  to  the  development  of  a technique  for  estimating  munition 
expenditure  requirements  for  the  specified  mission. 

Since  the  half-width  function  H(t)  was  defined  earlier  to  be  independent  of  the  wind 
direction/Fne-of-sight  angle,  the  reference  system  for  this  methodology  can  be  defined  to  facilitate 
reader  understanding.  The  front  to  be  screened  is  defined  parallel  to  the  x-axis  with  x = 0 being  the 
upwind  edge  of  the  front.  The  positive  y-axis  is  normal  to  the  front  in  the  upwind  direction.  The 
components  of  the  wind  vector  are  then  expressible  in  terms  of  their  relation  to  the  front;  i.e.. 


= u cos  6 

(2) 

= u sin  9 

(3) 

This  reference  system  provides  the  ability  to  identify  aiming  points  as  being  upwind  or  downwind 
of  the  front.  The  distance  between  the  two  fronts  is  the  range  (R). 

The  formation  of  a smoke  screen  is  a dynamic  process.  The  mission  requires  the  front 
to  be  screened  within  t0  seconds.  If  uy  > 0,  the  screen  must  be  formed  before  the  clouds  travel 
halfway  to  the  front  to  allow  time  for  replenishment  without  gaps  in  the  screen.  Since  the  effective 
width  of  a cloud  reaches  a maximum  at  some  time  tm,  the  screen  must  be  formed  prior  to  that 
time.  Hence  the  maximum  allowable  time  t*  to  form  the  screen  is  the  least  of  the  three  times;  i.e.. 


t*  = min 


f R 

jl0’  V “ 


100 


2u, 


(4) 


The  minimum  number  of  rounds  required  to  establish  the  screen  can  now  be  estimated  from  the 
width  of  a single  cloud  at  time  t*,  which  is  2H(t*).  The  minimum  number  is 


N = 


FRONT 

2H(t*) 


(5) 


which  must  be  increased  by  one  if  there  is  a fractional  round  required.  At  this  point,  one  must 
consider  the  constraints.  If  N > NTUB,  more  than  one  salvo  must  be  fired  to  establish.  For  the 
moment,  let  us  proceed  with  the  N < NTUB  case.  We  know  that  a single  salvo  will  establish  the 
screen.  This  is  the  simplest  case  for  which  the  replenishment  of  the  screen  can  be  considered.  The 
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initial  set  of  aiming  point  coordinates  for  NSE  establishment  salvos  are  expressed  as 
(Xj,yj):i=l,2, . . N,  where 


f 


I 


J 


Xj=  FRONT  - (2i  - 1)  [H(t*)J  -uxt* 


(6) 


yj  = uy  [t*  + TI  (NSE- 1 )] 

The  yj  values  are  the  minimum  setback  distances  required  to  place  the  screen  between  the  fronts. 
The  time  TI  is  the  interval  between  salvos  when  NSE  > 1.  In  order  to  maintain  the  screen  for  ST 
seconds,  it  is  necessary  to  fire  NTUB'  rounds  every  TI  seconds.  If  tm  is  sufficiently  large,  NTUB' 
may  be  less  than  NTUB,  thus  freeing  some  delivery  systems  from  the  smoke  mission.  The  trade-off 
between  choosing  TI  large  or  NTUB'  small  is  difficult.  Since  at  t = t*  the  most  upwind  cloud  is 
passing  the  upwind  edge  of  the  front,  TI  < t*  is  essential.  Let  us  choose  TI  = t*/NSR  so  that  as  few 
as  N/NSR  clouds  must  be  replaced.  Since  NSR  = 2 yields  the  largest  TI  = t*/2  and 
NTUB'  = ( 1 + sin  0)  N/2,  the  only  remaining  question  is  whether  the  cloud  starting  at  will  screen 
the  downwind  edge  of  the  front  when  it  gets  there.  This  is  assured  when  the  minimum  setback 
distance  y = Uy  (t*  + NSR  • TI)  and 


FRONT  - xN 


< t 


m 


(7) 


If  equation  7 is  not  satisfied,  then  the  most  upwind  aiming  point  that  satisfies  equation  7 must  be 
found  by  finding  the  xj  so  that 


FRONT -(2i-l)H(t*)-uxt* 


< t 


m 


(8) 


It  must  be  noted  also  that  NSR  can  be  chosen  as  2 only  if  t*/TMU  > 2 since  the  unit’s  capability  is 
a limiting  factor.  In  some  cases,  therefore,  NSR  = 1 is  required.  If  all  necessary  conditions  are 
satisfied,  then 


NTUB'  = max 


, x , . , [ 1 + NSR  - 1)  sin  0) 

1,  N - i + 

NSR 


The  maximum  duration  of  the  screen  is  determined  from  the  constraints  to  be 


T 


max 


= min 


j t*  LD-N  TI 
1 TMU  ’ NTUB'  ‘ 


(9) 


(10) 


12 


Li 

. 

f 

• ! 


I 


Now  let  us  turn  to  the  case  of  NSE  > I;  i.e.,  NTUB  < N.  Of  primary  concern  is  whether  the  fire 
unit  is  capable  ot  accomplishing  the  smoke  mission.  The  minimum  time  to  establish  is  t*  just  as  in 
the  NSE  = 1 case.  At  the  maximum  fire  rate,  the  maximum  possible  value  of  NSE  is 

NSEmax  = 11  1 + (II) 

Determination  of  the  required  value  of  NSE  must  consider  the  staggered  arrival  time  of  the  salvos 
and  the  possibility  that  some  aiming  points  may  not  require  repeated  firing.  With  the  maximum  fire- 
rate  assumed,  the  time  between  salvos  is  TMU.  In  this  interval,  a cloud  travels  a distance  of 
ux  • TMU.  The  separation  between  aiming  points  is  2H[t*  - TMU  • (NSE-1)],  which  arises  from  the 
fact  that  the  last  clouds  to  arrive  have  only  t*  - TMU  • (NSE-1)  seconds  to  grow  before  reaching 
the  upwind  edge  of  the  front.  The  aiming  points  are  held  constant  to  permit  the  maximum  fire  rate, 
if  required.  Denote  by  h the  value  of  H[t*-TMU  • (NSE-  1)J . Then  each  salvo  after  the  first  could 
produce 


NOLAP  = 


ux  • TMU 
2h 


(12) 


duplicate  clouds.  The  number  of  salvos  required  to  establish  the  screen  is  determined  by  solving 


NSE 


FRONT 

2 [(NTUB  - NOLAP)  • NSE  + NOLAP]  h 


The  equivalent  problem  to  solve  is  G(NSE)  = 0,  where 


(13) 


G(NSE)  - [(NTUB  - NOLAP)NSE  + NOLAP]  H [t*  - TMU  ■ (NSE  -1)]  - .5  FRONT  (14) 

If  G(NSEmax)  < 0,  then  the  unit  cannot  accomplish  the  mission.  If  G(NSEmax)  = 0, 

NSE  = NSEmax.  If  G(NSEmax)  > 0,  then  the  smallest  value  of  NSE  for  which  G(NSE)  > 0 should 
be  used. 


If  the  inequality 


H[t*-  TMU  -(NSE-1)] 


FRONT 

2[(NTUB  - NOLAP)  NSE  + NOLAP] 


> 0 


(15) 


is  not  satisfied,  the  unit  cannot  accomplish  the  mission.  Also,  if  NOLAP  > NTUB,  the  mission 
cannot  be  accomplished  because  the  unit  can  never  replace  all  the  clouds  required. 


It  should  be  recognized  that  if  NSE  < NSEmax,  then  the  maximum  fire  rate  may  not 
be  required  for  the  mission.  The  time  between  salvos  could  be  increased  to 


T1  = max[TMU,  t*/NSE] 


and  the  maximum  time  this  rate  of  fire  could  be  maintained  would  be  increased  theoretically  to 


TUM'  = max 


based  on  a constant  number  of  rounds  fired. 

A consideration  of  the  maintenance  of  the  screen  takes  the  above  results  into  account. 
At  the  maximum  fire  rate,  the  number  of  replenishment  salvos  available  is  (NSEmax  - NSE), 
denoted  by  NSRmax.  The  adjustment  to  TUM  may  be  sufficient  to  increase  NSEmax  and  is  made 
by  substituting  TUM'  for  TUM  in  equation  1 1 . The  new  value  of  NSRmax  is 


NSRmax  NSRmax  + 


TUM’ 


A further  modification  of  NSRmax  is  possible  if  NTUB  > NTUB'.  Then  one  can  assume  that  more 
than  one  delivery  system  (tube)  can  be  aimed  at  a particular  aiming  point  and  fired  in  rotation  to 
replenish  the  screen.  Hence,  NSRmax  can  be  adjusted  by  a factor  of  NTUB/NTUB'.  The  maximum 
duration  of  the  screen  is  obtained  by  considering  the  limiting  factors  of  performance  and  logistics  in 
the  expression 


Tmax  = TI  • min  )NSRmax 


NTUB 

NTUB' 


LD  - NSE  • NTUB 
NTUB' 


Several  cases  were  noted  above  in  which  the  fire  unit  is  not  capable  of  establishing  the 
required  screen  in  the  t*  seconds  allowed.  Rather  than  dismissing  the  mission  as  impossible,  an 
analyst  would  prefer  to  examine  the  problem  with  a weaker  constraint  on  the  value  of  t*,  which  is 
usually  specified  arbitrarily.  The  only  constraint  that  will  be  imposed  on  t*  is  that  t*  < tm.  If 
t*  > tm,  the  clouds  would  be  decreasing  in  width  at  the  time  they  must  be  merging.  The  minimum 
possible  value  of  t*  must  satisfy  equation  15.  Of  course,  if  t*  = tm  does  not  satisfy  equation  15, 
then  it  is  reasonable  to  conclude  the  unit  cannot  provide  the  required  screen.  If  a t*  < tm  does 
exist,  then  NSEmax  becomes 


NSEmax  = 1 + 


min  (TUM,  tm  - t*) 


and  NSR|)iax  cun  be  adjusted  accordingly.  If  NSK  > NSEmax  even  after  this  adjustment,  then  the 
unit  cannot  provide  the  screen.  Otherwise,  the  procedure  followed  with  the  initially  specified  t*  can 
be  used  to  estimate  the  expenditure  requirements. 

Once  the  values  of  NSE,  NTUB',  NOLAP,  and  t*  have  been  determined  by  the 
appropriate  method,  it  is  necessary  to  determine  the  number  of  rounds  required  both  to  establish 
the  screen  and  to  maintain  it  for  ST  seconds.  It  is  also  possible  to  specify  a set  of  aiming  points 
which  will  accomplish  the  mission.  For  NSE  > 1,  NTUB’  = NTUB,  but  not  all  tubes  repeat  fire  at  a 
single  aiming  point.  Some  tubes  must  fire  at  new  aiming  points  determined  by  the  value  of  NOLAP. 
Consequently,  we  must  specify  the  coordinates  for  the  i-th  salvo  and  the  j-th  tube.  The  Xy  values  are 
given  by 


U 


x j j = FRONT  - (2  k - 1)  H [t*  - TI(NSE)]  - ux  t* 


where 


k = j • NTUB  - (j  * 1)  NOLAP 


The  corresponding  yy  values  are 


yy  = uy  [t*-TI  (NSE  - 1)] 


If  NSE  = 1,  only  the  yy  values  are  different,  being  given  by 


Vy  = 2 Uy  t* 


but  the  aiming  points  for  replenishment  fire  are  only  those  for  which 

j > NTUB  - NTUB' 


(21) 


(22) 


(23) 


(24) 


(25) 


If  the  required  duration  ST  of  the  screen  is  less  than  the  maximum  possible,  Tmax,  the  number  of 
rounds  required  to  accomplish  the  mission  is 


CT 

Nt  = NRE  + NTUB'  • II  — 
1 TI 


(26) 


where  the  number  of  rounds  to  establish  is 


NRE  = NSE  • NTUB 


(27) 


This  methodology  for  instantaneous  sources  has  been  implemented  in  a Fortran  V code 
called  MUNEXP,  which  is  described  elsewhere  in  this  report. 


B.  Continuous  Sources. 


1 . Effective  Screen  Dimensions. 

The  lack  of  polar  symmetry  in  continuous  source  plumes  requires  separate 
methodologies  for  the  two  primary  line-of-sight  orientations. 

The  only  successful  method  to  date  for  estimating  crosswind  screen  lengths  is  a 
systematic  search  procedure  based  upon  experience  in  prior  efforts.  If  the  emission  rate  is  constant, 
which  is  assumed  frequently  in  the  absence  of  emission  rate  measurements,  the  screen  length  grows 
with  time  to  a maximum  length  and  reaches  a steady  state  until  emission  stops.  After  burnout  at 
t = tg,  the  upwind  end  of  the  screen  moves  downwind  until  the  screen  vanishes.  If  the  emission  rate 
is  not  constant,  the  cloud  must  be  simulated  for  an  array  of  times  to  get  the  time  trace  of  the  ends 
of  the  screen.  It  is  likely  that  the  screen  length  will  stabilize  at  some  value  during  the  burning  time. 
If  multiple  sources  are  involved,  stabilization  is  more  likely  since  the  cumulative  effect  of  many 
clouds  at  different  stages  of  growth  tends  to  smooth  the  variations  of  individual  clouds.  The 
recommended  procedure  is  to  determine  the  screen  length  just  prior  to  burnout.  In  most  cases,  the 
result  is  a screen  length  that  is  relatively  stable  and  has  existed  for  a large  portion  of  the  burning 
time.  This  screen  length  is  the  most  significant  measure  of  effectiveness  for  a munition  because  it 
usually  forms  in  a short  time  and  has  significant  durability  and  its  location  relative  to  the  center  of 
the  impact  zone  is  predictable  for  use  by  mission  planners.  The  formation  time  can  be  estimated  by 
dividing  the  distance  to  the  downwind  end  of  the  screen  by  the  average  windspeed. 

The  location  of  the  two  ends  of  the  screen  should  be  sought  independently.  By 
choosing  a small  set  of  x distances  which  span  a reasonable  interval  in  which  an  end  may  be 
expected  to  lie,  one  can  use  a small  number  of  simulations  combined  with  an  appropriate 
interpolation  scheme  to  find  the  end  sought.  The  expense  of  this  procedure  depends  strongly  on  the 
accuracy  required,  but  the  accuracy  of  the  simulation  model  should  be  considered  before  imposing 
resolution  criteria. 

For  a headwind/tailwind  situation,  the  effective  screen  width  of  a single  continuous 
source  is  much  smaller  than  in  the  crosswind  case.  If  the  screen  must  be  established  within  a 
specified  minimum  time,  the  width  at  that  time  will  likely  be  the  minimum  value  for  all  later  times 
less  than  the  time  at  which  the  smoke  emitted  at  t = tg  reaches  the  front  to  be  screened.  Estimation 
of  the  effective  screen  width  is  accomplished  by  using  the  instantaneous  source  technique  described 
earlier. 

2.  Munition  Expenditure  Requirements. 
a.  Crosswind  Line  of  Sight. 

Once  an  effective  crosswind  screen  length  is  determined  for  a single  munition  or 
perhaps  several  munitions  fired  at  the  same  aiming  point,  the  construction  of  a screen  to  accomplish 
the  mission  must  be  undertaken.  Again,  because  munition  effectiveness  is  the  prime  objective,  it  is 
appropriate  that  the  ability  of  troops  to  hit  an  aiming  point  be  ignored.  Therefore,  a smoke  screen 
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can  theoretically  be  constructed  by  placing  smaller  screen  segments  end  to  end.  To  avoid  possible 
gaps  in  the  screen,  only  90  percent  of  the  length  of  an  effective  segment  is  used,  implying  overlap  of 
segments.  Thus,  if  a screen  of  length  FRONT  is  to  be  built  of  segments  having  length  S,  the 
minimum  number  of  segments  (or  aiming  points)  required  is 

No  =1  + 11  — II  (28) 

3 ,9S 

When  a single  round  produces  the  screen  segment  S,  the  number  of  rounds  required  equals  the 
number  of  segments.  However,  some  small  munitions,  such  as  mortars,  do  not  produce  significant 
screens  with  a single  round.  In  these  cases,  the  segment  length  ought  to  be  estimated  for  multiple 
rounds  fired  at  a single  aiming  point.  Delivery  errors  should  be  a part  of  this  type  of  simulation. 
Operational  considerations  should  drive  the  trade-off  between  the  number  of  aiming  points  and  the 
total  number  of  rounds  required  toward  the  selection  of  fewer  aiming  points. 

Long  burning  times  of  continuous  sources  usually  do  not  cause  conflict  with  the  fire 
rate  of  the  unit.  Therefore,  munition  expenditure  requirements  are  obtained  as  the  product  of  Ng 
and  the  number  of  rounds  per  segment,  Ng/g: 

NRE  = Ng  • N j^/g  ( 29 ) 

Maintenance  of  the  screen  is  achieved  by  repeating  the  establishment  fire  at  a time 
interval  based  on  the  burn  time  and  formation  time.  Suppose  the  upwind  end  of  the  screen  is 
formed  at  tf  seconds  after  function  time.  Then  the  screen  exists  for  tg  - tf  seconds,  after  which  the 
edge  begins  to  move  downwind.  If  tg  seconds  are  allowed  between  fires,  there  is  a possibility  of  a 
gap  in  the  screen.  If  tf  is  small,  say  tf  < 10  seconds,  the  combination  of  the  10  percent  overlap  and 
the  arrival  of  clouds  from  upwind  sources  will  probably  suffice  to  allow  the  selection  of  tg  as  the 
replenishment  interval.  If  tf  cannot  be  ignored,  let  tg  - tf  be  the  interval.  The  analyst  may  have  to 
improvise  if  the  screen  length  is  not  stable  due  to  variable  emission  rates. 

If  the  screen  is  to  be  maintained  for  a period  of  ST  seconds,  the  total  number  of 
rounds  required  is 


Nt  = NRE 


+||_ST_ 

\ tg  - tf 


(30) 


This  methodology  is  the  simplest  in  concept  and  application;  hence,  no  software  is  included  in  this 
report. 


b.  Headwind/Tailwind  Line  of  Sight. 

The  problem  of  creating  a screen  under  headwind/tailwind  conditions  is  the  most 
difficult.  Fortunately,  the  frequency  of  occurrence  in  the  real  world  of  this  situation  is  very  low. 
Since  each  source  produces  a short  screen  at  some  specified  time,  the  lateral  distance  between 
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sources  is  critical  to  the  avoidance  of  clear  “lanes”  between  plumes.  For  smoke  generators  with 
operators,  it  is  a simple  matter  to  place  multiple  generators  at  ideal  sites  for  the  desired  screening 
effect.  The  methodology  described  here  considers  only  rounds  with  multiple  submunitions  with 
random  distribution  of  the  impact  points. 


Assume  the  half-width  of  the  screen  produced  by  a single  submunition  is  w, 
determined  by  the  method  in  section  III,A,1.  Let  x represent  the  crosswind  axis  with  the  origin 
being  the  center  of  the  front  to  be  screened.  If  all  rounds  are  fired  at  the  point  (0,y),  where  y is  the 
upwind  distance  from  the  front  required  to  place  the  screen  between  the  target  and  observer,  the 
only  distribution  of  impact  points  to  be  considered  is  that  with  respect  to  the  x-axis. 

For  each  round,  we  will  assume  n submunitions  are  distributed  uniformly  in  an  area 
bounded  by  a circle  of  radius  R.  The  distribution  of  impact  points  with  respect  to  x is  given  by  the 
probability  density  function  (pdf) 


fix)  = 


2 \/  R2  - x2  for  lxl<  R 


for  I x I > R 


which  assumes  the  center  of  the  impact  area  to  be  at  x = 0.  It  is  preferable  to  rewrite  this  equation 
to  represent  the  distribution  of  impact  points  u (parallel  to  x)  centered  at  v: 


f(u  - v)  = 


i — ~ [R2-(u-v)2]  2 for  I u - v I < R 
IttR2 


for  I u - v I > R 


Let  the  front  be  represented  as  a set  of  points  x lying  in  the  interval  (-  FRONT/2,  FRONT/2).  The 
probability  that  an  arbitrary  point  x is  screened  by  one  of  the  submunitions  distributed  by 
equation  32  is  the  probability  that  the  distance  between  x and  impact  point  u is  less  than  the  screen 
width  w;  i.e., 

Pj  = P [ 1 x - u I < w] 

= P [ lu  - xl<  w] 

= P [-w  < u - x < w] 

(33) 

= P[x-w<u<x  + w] 


/*X  + ' 

J\-\ 


f (u  - v)  du 


Lot  y = u - v to  pet 


!X  T 

X - 


x + w - v 


f(y)dy 


X - W - V 


If  equation  31  is  used  with  equation  34,  the  result  after  integration  is 


?!  (x,  v)  = G (x  + w - v)  - G (x  - w - v) 


where 


G(x)  = 


\ f(y)  dy  = — [ 
ttR2  l 


x /r2  - x2  + R2  sin  1 


(*>]' 


if  I x I < R 


if  x > R 


if  x < -R 


The  probability  that  the  point  x is  screened  by  at  least  one  of  the  n submunitions  is 

Pn(x,v)  = 1 - [1  - Pj  (x,v)]n  (37) 

If  m rounds  are  fired  at  the  aiming  point,  the  probability  of  the  point  x being  screened  by  at  least 
one  of  the  mn  submunitions  is 


Pnm^x)  ^ Pn  ^x’ vi^ 


The  form  of  equation  38  is  too  involved  to  permit  an  analytic  evaluation  of  Pnm(x)  or  its  expected 
value,  which  represents  the  fraction  of  the  front  which  would  be  screened  by  m rounds.  A Monte 
Carlo  simulation  technique  can  be  used  to  evaluate 


E (FRONT, n,m,w)  = 


FRONT/2 


FRONT 


Pnm<x>  dx 
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E quation  40  is  the  variance  ot  the  expected  screening  coverage  li. 

The  method  used,  which  has  been  implemented  in  the  computer  program  MUNX 
(described  later),  simulates  repeated  firings  of  nt  munitions  at  the  aiming  point  and  averages  the 
results  ot  the  many  experiments.  The  value  of  m is  increased  until  E exceeds  a specified  value. 

C.  Comments  on  Methodology. 


All  of  the  techniques  described  above  have  the  characteristic  of  overestimating 
expenditure  requirements  to  an  unknown  degree  because  they  neglect  the  screening  effect  of 
overlapping  cloud  regions  outside  the  arbitrary  boundary.  This  is  done  intentionally  for  several 
reasons : 

1 . To  simplify  analysis. 

2.  To  compensate  for  assumptions  of  ideal  round  delivery. 

3.  To  prevent  unrealistic  comparisons  between  munition  systems  in  competitive 

studies. 


The  methodologies  are  based  on  idealized  scenarios.  Their  validity  in  representing  the 
real  world  has  not  beer,  examined.  It  is  assumed  that  their  application  would  be  valid  for 
comparative  effectiveness  studies  in  which  absolute  effectiveness  estimates  are  not  required. 

IV.  COMPUTER  PROGRAMS. 

A.  MUNEXP. 

The  program  MUNEXP  implements  the  methodology  of  section  III,  A,  for 
instantaneous  sources.  The  code  is  in  Univac  Fortran  V and  requires  8,357  words  of  user  memory'. 
The  program  is  constructed  in  modular  form  to  facilitate  user  modification  of  the  half-width  model. 
The  input  format  is  free-field  to  permit  easy  use  of  a demand  terminal.  For  the  same  purpose,  all 
output  is  restricted  to  70-character  lines. 

Consider  as  an  example  the  following  screening  mission.  A 200-meter  front  must  be 
screened  below  10  percent  transmittance  within  45  seconds  of  the  arrival  of  the  first  of  many 
1 55-mm  Ml  10  WP  smoke  projectiles.  The  wind  is  blowing  toward  the  front  at  an  angle  of  30°  and  a 
speed  of  10  mph.  The  meteorological  stability  is  Pasquill  category  C (neutral).  The  relative  humidity 
is  50  percent.  The  mission  is  to  be  accomplished  by  two  artillery  batteries  with  eight  tubes  each. 
Each  battery  carries  200  rounds  into  the  field.  The  crews  are  capable  of  firing  one  round  every 
15  seconds  for  up  to  4 minutes.  The  screen,  once  established,  should  be  maintained  for  4 to 
8 minutes.  The  estimated  target-observer  range  is  1 500  meters.  The  questions  to  be  addressed  by  the 
analyst  are : 

1 . Can  the  mission  be  accomplished? 

2.  If  so,  how  many  rounds  are  required? 
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3. 


Where  should  the  rounds  be  aimed? 


4.  How  ean  the  screen  be  maintained? 

5.  How  long  can  the  screen  be  maintained? 

The  first  step  in  the  analysis  is  the  estimation  of  the  half-width  of  a single  round.  The 
round  is  simulated  under  the  given  meteorological  conditions  to  generate  a matrix  of  CL  values, 
shown  in  table  1,  with  the  interpolated  estimates  of  the  half-width  corresponding  to  a CL  of 
693  mg/m-  shown  in  the  augmented  column.  For  WP,  this  CL  will  give  a 10  percent  transmittance 
in  the  visible  wavelengths. 

Using  the  x*  values  of  table  1 as  the  half-widths  at  the  times  given,  a least-squares  fit 
with  a 4-th  degree  polynomial  model  produces  the  curve  shown  in  figure  1.  From  this  curve,  one 
can  estimate  that  the  maximum  width  occurs  at  t = 120,  which  is  the  input  value  for  TM.  The 
equation  of  the  model  is  given  at  the  top  of  figure  1,  the  coefficients  being  the  A(l)  input. 

Applying  all  the  information  given  to  the  computer  program  MUNLXP  produces  the 
output  shown  in  figure  2,  which  includes  the  input  data  and  the  solution.  The  output  indicates 
three  rounds  are  to  be  fired  at  the  required  aiming  points,  but  only  the  third  point  is  required  for 
replenishment  every  22.5  seconds.  The  screen  can  be  maintained  for  90  minutes,  requiring  a total  of 
1 3 to  24  rounds  for  a 4-  to  8-minute  screen. 

Most  variable  names  in  MUNEXP  correspond  to  the  notation  of  this  report,  but  the 
code  may  not  always  use  a variable  name  exclusively  as  used  in  the  text.  Table  2 gives  the  input 
variables,  format,  and  definitions  required  for  using  the  program.  Table  3 shows  the  external 
subroutines  required  and  their  function.  The  listing  of  the  code  in  figure  3 includes  the  routines 
used  for  the  example  problem. 

B.  MUNX. 

The  program  MUNX  implements  the  methodology  described  in  section  III.B.2  for 
continuous  sources  with  a headwind/tailwind  line  of  sight.  The  code  is  written  in  Univac  Fortran  V 
and  requires  8800  words  of  user  memory.  Like  MUNEXP,  the  I/O  format  is  designed  for  user 
convenience  and  use  of  demand  terminals;  input  is  free  field  and  output  is  limited  to  70  characters. 

For  an  example  problem,  consider  the  requirements  to  form  a screen  125  meters  long 
using  the  1 55-mm  Ml  16  HC  projectile.  The  screen  must  be  formed  45  seconds  after  impact.  Under 
10  mph,  category  B (lapse)  conditions,  each  of  four  submunitions  produces  a 3.8-meter  screen  at 
45  seconds.  The  delivery  errors  are  assumed  to  be  arange  = 4 and  ^deflection  = 3.  The  four 
submunitions  fall  within  74  meters  of  the  impact  area  center  95  percent  of  the  time. 

To  get  a 125-meter  screen  with  100  percent  probability  is  impossible  under  this 
methodology.  It  is  reasonable  to  ask  for  95  percent  certainty  and  not  unreasonable  to  think  the 
corresponding  real  world  experiment  would  give  a larger  screen.  Hence,  we  will  choose  EMIN  = .95. 
We  will  seek  a precision  level  of  three  significant  digits  in  the  integration  procedure.  This  requires 
EPS  to  be  5 X 10"4.  To  prevent  excessive  run  time,  we  will  choose  NREP  = 30,  NMAX  = 20,  and 
MMAX  = 30. 


2! 


Table  1.  Matrix  of  CL  Values 


5ET  1.  0-2-59PCT 

V=  ^5.54  * «»  1 -c..  liil-.Jl-t  + •‘-2e-05X**3 


Half-Width  Curve 


[ 

f 


• 5XQT  F.MEXP 


WINDSPEED.  M/SEC 
NINO  ANGLE,  DEG 
RANGE.  M 
FRONT,  M 


NUMBER  OF  TUBES  AVAILABLE 
NUMBER  OF  ROUNDS  AVAILABLE 
MINIMUM  TIME  BETWEEN  SALVOS,  SEC 
MAX  DURATION  OF  MAX  FIRE  RATE,  SEC 


4.470 

30.00 

1500.00 

200.00 


16 

400 

15.00 

240.00 


SPECIFIED  TIME  TO  ESTABLISH,  SEC  45.00 
PREDICTED  TIME  TO  ESTABLISH,  SEC  45.00 
ROUNDS  TO  ESTABLISH  3 
ROUNDS  FOR  REPLENISHMENT  1 
TIME  TO  REPLENISH,  SEC  22.50 


MAXIMUM  SCREEN'DURATION,  MIN 


90.00 


REQUIRED  AIMING  POINTS 


INDEX 

X Y ■ MIN 

SETBACK 

REPLENISH(o) 

1 

-1B.B35 

201 .150 

2 

-108.104 

201.150 

3 

-197.372 

201 .150 

• 

i 


ROUNDS 

REQUIRED 

FOR 

4.0-MIN 

SCREEN 

13 

ROUNDS 

REQUIRED 

FOR 

6.0-MIN 

SCREEN 

18 

ROUNDS 

REQUIRED 

FOR 

8.0-MIN 

SCREEN 

24 

Figure  2.  MUNEXP  Program  Output 


Table  2.  MUNHXP  Input 


Card  No. 

Variable 

Type 

Def lnition 

Integer 

Integer 


Integer 

Integer 

Real 

Real 


Number  of  parameters  in  half-width  model. 
Number  of  screen  duration  times. 

Parameters  of  half-width  model. 

Estimate  of  time  of  maximum  half-width,  sec 

Length  of  front,  m. 

Wlndspeed,  m/sec. 

Required  minimum  formation  time,  sec. 

Wind  angle  with  front,  degree  < 90°. 
Observer-target  distance,  m. 

Number  of  available  delivery  systems; 
default  value  is  1000. 

Number  of  available  rounds;  default  value 
is  10000. 

Minimum  firing  interval,  sec;  default 
value  is  1. 

Maximum  duration  of  maximum  fire  rate, 
sec;  default  value  is  1000. 

Screen  duration  times,  min. 


J=1 ,2 , . . KT 

- 

‘Cases  can  be  stacked  by  repeating  cards  3-5. 


Table  3.  MUNEXP  Subroutines 


Name 

Call 

Function 

HW 

F=HW(T) 

Computes  value  of  half-width  model 
at  time  T. 

DHW 

F=DHW(T) 

Computes  derivative  of  half-width 
function  at  time  T. 

HTI 

F=HTI(T) 

Evaluates  left-hand-side  of  equation  15. 

NR 

CALL  NR(X,FCT) 

Newton- Raphson  routine  to  solve 

FCT(X)  = 0. 

' 

1 1 


DRDIFF 


CALL  DRDIFF 
(DYDX,  FCT,  X) 


Evaluates  derivative  of  FCT  at  X. 


SAD*  CMOKEX ( 1 ) 
1 
2 

3 

4 

5 

6 
7 

a 

9 

10 

11 

12 

13 

14 

15 

16 
17 
16 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 


MUNEXP 

D I MENS  I ON  ST ( 10) ,X(200), 7(200) 

COMMON /ONE /M , A ( 1 0 ) 

COMMON/TWO/  IN, FRONT 

DATA  NTUBE , LOD, T TMU . TTUM/ 1 000 , 1 0000 , 1 . , 1 000 . / 

EXTERNAL  HW.DHW.HTI 

READ! 5,100)  M,KT 

RE AO ( 5 , 1 00 ) (A(I),  1=1 , M ! , i M 

10  READ(5,100,END=99)  FPONT , U , TO , DANG , RNGE 
RE AD ( 5 , 100)  NTUB.LD.TMU.TUM 

IF(NTUB. LE.O)  NTUB=NTU8E 
IF(LD.LE.O)  LD=LOD 
IF(TMU. LE.O.  ) TMU=  TTMU 
IF(TUM. LE.O.  ) TUM=  TTUM 
WR ITE(6, 105)  U, DANG, RNGE, FRONT 
WR I TE( 6 , 1 1 1 ) NTUB.LD.TMU.TUM 
RANG=DANG*3. 1415927/180. 

' UX=U*COS(RANG) 

' UY=U*SIN( RANG) 

ST ( 1 ) = 1 0 . 

IF(KT.GT.O)  READ! 5,100)  ( ST ( I ) , 1=1, KT) 
TP=AMIN1 (TO, (RNGE- 100. )/(2*UY+.01 ) ) 

CALL  NR(TM.DHW) 

TSTAR=AMIN1(TP,TM) 

AN= . 5* FRONT/HW( TST  AR ) 

N = AN 

I F ( AN-N.GT.O.  ) N = N + 1 
METH=2 

NSEMAX= 1+TST  AR/TMU 
IF(N.LE.NTUB)  ME  TH  = 1 
IF(METH.EQ. 1 ) GO  TO  50 
AN  = N 

NSE= AN/NTUB+ . 999 
NSET  =NSEMAX+ 1 

11  NSET  =NSET- 1 

HT  = HW( TST  A R- TMU* ( N SET- 1 ) ) 

NOLA  P= . 5*UX» TMU/HT 
NP=NSET* (NTUB-NOLAP)+NOLAP 
G=(NOLAP+( NT  UB-NOL  AP ) *NSET ) *HT" • 5*  FRONT 
I F ( G I 13,12,11 

13  IF(NSET.GE.NSEMAX)  GO  TO  30 
NSET  =NSET  + 1 

12  HT=HW( TSTAR- TMU* (NSET— 1 ) ) 

NOLA  P= . 5* UX*  TMU/HT 
NP=NSET*(NTUB-NOLAP)+NOLAP 

I F ( HT . LT. .5*  FRONT/NP ) GO  TO  30 
TP  = 0 

NSE=NSET 

14  NIMAX=NSEMAX-NSE 
TS=TMU* ( NSE- 1 )+TP 
TSTAR=AMAX1 ( TS, TO) 

TI=AMAX1  (TMIJ.TSTAR/NSE) 

TUMP  = TUM*  T I / TMU 

NSEP= 1 +M IN ( TUMP , TM-TP ) /TMU 
NIMAX=N!MAX+(NSEP-NSEMAX) 
HT=HW(TSTAR-TI*(NSE-1 ) ) 

NTUBP=NTUB 

IF(NOLAP.GE.NTUB)  GO  TO  40 


Figure  3.  MUNEXP  Program  Listing 


I 


t . 


t 


H 


L 

J 


GO  TO  70 
30  TP=TM*.9 
IN  = NP 

WRITEC6, 130) 

IF(HTI (TM) . LT.O. ) GO  TO  40 
CALL  NR(TP.HTI) 

TP=TP+TMU* (NSE-1 ) 

IF(TP.GT.TM)  GO  TO  40 
NSEMAX.1+AMIN1 ( TUM , TM-TP) /TMU 
T P = T P-TMU* (NSE-1  ) 

NSE*NSET 

IF(NSE.LE.NSEMAX)  GO  TO  14 
40  WR ITE( 6 , 140) 

GO  TO  10 

50  NIMAX*TUM/TMU-1 
NSE=1 

NSR=  1 

I F ( N IMAX . GT . 1 ) NSR  = 2 
HT=HW( TSTAR) 

1 = AN 

51  TP»(  FRONT- (2*  1-1  ) *HT—UX* TSTAR) / ( UX+-1  ,E-5) 
IF(TP.LE.TM)  GO  TO  52 

1 = 1-1 

I F ( I . LE . 1 ) GO  TO  52 
GO  TO  51 

52  NTUBP=N-I+( 1+fNSR-l ) *S IN( RANG) ) /NSR 
NTUBP=MAX ( 1 , NTUBP) 

T I = T STAR/NSR 

NIMAX=NIMAX*NTUB/NTUBP 

NTUB=N 

70  NIMAX=MIN(NIMAX, (LO-NSe*NTUB)/NTUBP) 
TMAX=NIMAX*T 1/60 . 

NRE=NSE*NTUB 

WRITE (6, 104)  TO. TSTAR, NRE, NTUBP, TI 
WR I T E( 6 , 1 1 0 ) TMAX 
I F ( NSE . EO • 1 ) GO  TO  80 
WR I TE( 6 , 1 1 4 ) 

GO  TO  81 

80  WRITEC6.107) 

81  K = 0 

IF(NOLAP.EO.NTUB)  NSE*1 


99 

DO  86  L= 1 , NSE 

100 

DO  85  d = 1 , NT UB 

101 

K = K+1 

102 

REP= ' 1 

103 

IF(J.GT.NTUB-NTUBP)  REP='*' 

104 

X( J)=FRONT-( 2*K-1 ) *HT-UX* TSTAR 

105 

Y ( d ) =UY* ( T ST  AR+T I * ( NSE— 1 )) 

106 

I F ( METH . EQ . 1 ) Y(d)=2*UY*TSTAR 

107 

IF(NSE.GT. 1 ) GO  TO  83 

108 

WR I TE( 6 , 1 08 ) d,X(d) ,Y(J) ,REP 

105 

GO  TO  85 

110 

83 

WR I TE ( 6 ■ 1 1 3 ) L,d,X(d),Y(d),REP 

111 

85 

CONTINUE 

1 12 

K=K-NOLAP 

1 13 

86 

CONT INUE 

114 

WRITEC 6,112) 

115 

DO  90  J* 1 , KT 

Figure  3.  (Continued) 
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116 
117 
1 IB 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 
1 38 
139 
1 40 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 


IF(ST( J(.GT.TMAX)  GO  TO  95 
K*60.*ST(J)/TI-,01 
NT=NRE+K*NTUBP 
WRITE16.106)  ST ( J ) , NT 
GO  TO  90 

95  WR I T E ( 6 , 109)  ST(J) 

90  CONTINUE 
GO  TO  10 
99  STOP 
100  FORMAT() 

105  F0RMAT(1H0,////,2X'WIN0  SPEED.  M/SEC ' . T60 , F 1 0 . 3 , / , 

* 2X 1 W I ND  ANGLE,  DEG  . ' , T60 , FI  0 . 2 , / , 

« 2X 1 RANGE , M'  .T60.F10.2,/, 

* 2X 1 FRONT , M'  , T60 , F 1 0 .2 , /// ) 

104  FORMAT ( 

*2X'SPECIFIED  TIME  TO  ESTABLISH,  SEC ' , T60 , FI  0 . 2 , / , 
*2X'PREDICTED  TIME  TO  ESTABLISH,  SEC ' , T60 , FI  0 .2 , /, 

* 2X  ' ROUNDS  TO  EST AB LI SH ' , T60 , 1 1 0 , / , 

* 2X ' ROUNDS  FOR  REP LENI SHMENT ' , T60 , 1 1 0 , / , 

•2X'TIME  TO  REPLENISH, SEC' .T60.F10, 2,///) 

106  FORMAT ( 2X 1 ROUNDS  REQUIRED  FOR'.FS.I,'  MIN.  SCREEN T60 , 1 1 0 ) 

107  FORMAT (////, T24 , 'REQUIRED  AIMING  POINTS',//, 

»2X' INDEX' , T1 7, 'X' . T24, ' Y*MIN  SETBACK ’ ,T47, ’ REPLENISH! * ) ' ,//) 

108  FORMAT (3X, 1 3 , 6X , FI  0 . 3 , 5X , FI 0 . 3 , 1 0X , A 1 ) 

109  F0RMAT(2X' ROUNDS  REQUIRED  FOR'.FS.I,'  MIN.  SCREEN', T60, 

* 'NOT  AVAILABLE' ) 

110  F0RMAT(2X'MAXIMUM  SCREEN  DURATION,  MIN.' , T60 , F 1 0 . 2 ) 

111  FORMAT! 

* 2X ' NUMBER  OF  TUBES  AVA I L ABLE ' , T60 , 1 1 0 , / 

* 2X ' NUMBER  OF  ROUNDS  AV A I LABLE ' , T60 , 1 1 0 , / 

* 2X ' M l N I MUM  TIME  BETWEEN  SALVOS,  SEC . ' , T60 , F 1 0 . 2 , / 

* 2X ' MAX  DURATION  OF  MAX  FIRE  RATE,  SEC . ' , T60 , F 1 0 . 2 , /// ) 

112  FORMAT!///) 

113  FORMAT! IB.8X, I3,6X,F10.3,5X,F10.3, 1 OX  , A 1 ) 

114  FORMAT!////, T24, 'REQUIRED  AIMING  POINTS',//, 

* 2X ' SALVO ' 5X ' INDEX'  ,T31  , 'X' ,T40,  ' Y*MI N SETBACK'  ,T58, 

* 1 REPLENISH! * ) ' ,//) 

130  F0RMAT(10X'SCREEN  CANNOT  BE  ESTABLISHED  BY  FIRE  UNIT  IN' 

* , ' SPECIFIED  TIME. ' ,//) 

140  FORMAT! lOX'SCREEN  EXCEEOS  CAPACITY  OF  FIRE  UNIT.') 

END 


Figure  3.  (Continued) 
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SAD* SMDKEX ( 1 ) 

.NEWRAP 

1 

SUBROUTINE  NR(X.FCT) 

2 

EXTERNAL  FCT 

3 

COMMON/C/ 1 ERR, F, N , NMAX 

4 

NMAX  = 30 

5 

I ERR  =0 

6 

XO  = X 

7 

EPS=5. E-3 

8 

N=  1 

9 

1 CALL  DRDIFF(DFDX.FCT.XO) 

10 

F * FCT ( XO ) 

11 

EMIN=1 ,E-4 

12 

I F ( ABS ( DFDX ) . LT, EMIN)  DFDX«S IGN ( EMIN , DFDX) 

13 

DX=-F/DFDX 

14 

X1.X0+DX 

15 

ERR=ABS(DX)/(ABS(X1 )+1 ,E-6) 

16 

IF(ERR.LE.EPS)  GO  TO  5 

17 

X0=X1 

18 

N = N+ 1 

19 

IF(N. LT.NMAX+1 ) GO  TO  1 

20 

F=FCT ( XI ) 

21 

IF(ABS(F).LT..1)  GO  TO  5 

22 

I ERR  = 1 

23 

N=N"  1 

24 

5 X*X1 

25 

RETURN 

26 

END 

SAD*SM0KEX(1 ).drdiff 

1 SUBROUTINE  DRDI FF( DFDX , FUNC , X) 

2 COMMON/D/HH 

3 DATA  NB/4/ 

4 C NB  = 1/2  THE  NUMBER  OF  DIGITS  CARRIED  BY  THE  MACHINE  BEING  USED 

5 F*  FUNC ( X ) 

6 XX*X+HH 

7 D=.001*ABS(XX) 

8 IF(ABS(XX) . LT. 10,**(3“NB) )D«10. »*-NB 

9 H*  AM  INI  ( F , D) 

10  DFDX* ( FUNC ( XX+H) “F ) /H 

11  HH=0 

12  RETURN 

13  END 


Figure  3.  (Continued) 


29 


SAD*SM0KEX(1  ) .WIDTH 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


1 


FUNCTION  HW( T) 

COMMON/ONE/M. A( 10) 

SUM= A( 1 ) 

DO  1 I =2  , M 

SUM=SUM+A( I ) *T** ( 1-1  } 
HW-SUM 
GO  TO  10 
ENTRY  DHW(T) 

SUM«A(2) 

DO  2 1=3  M 

2  SUM=SUM+A( I)*(I-1)*T**(I-2) 
HW*SUM 
10  RETURN 
END 


SAD*SMOKEX(  1 ) . FUNC  EX 

1 FUNCTION  HTI (T) 

2 CCMMON/TWO/'  IN,  FRONT 

3 HT I =HW( T )“ . 5*  FRONT /IN 

4 10  RETURN 

5 END 


Figure  3.  (Continued) 


Applying  this  information  to  the  program  MUNX  produces  the  output  in  figure  4, 
which  indicates  28  rounds  are  required  to  cover  more  than  95  percent  of  the  125-meter  front.  The 
variance  in  the  expected  coverage  is  7.54  X 10-4,  which  can  be  used  to  express  the  expected 
coverage  as  .953  ± .0275.  The  output  also  says  that  16  points  were  used  in  each  integration  to  get 
an  average  relative  error  of  6.94  X 10"5  in  the  coverage  estimate.  The  variance  within  the 
30  replications  is  shown  in  the  last  column  to  provide  insight  to  the  range  of  coverage  values 
observed. 


This  computer  run  required  nearly  3 minutes  of  CPU  time  on  the  Univacll08 
computer,  indicating  that  user  input  causing  premature  termination  could  be  expensive.  There  are 
ways  to  accelerate  the  procedure  which  will  be  implemented  in  future  work. 

Input  format  and  variables  are  shown  in  table  4.  Table  5 gives  the  external  subroutines 
and  their  functions.  Some  additional  comments  are  required  for  clarification.  The  methodology  is 
an  iterative  technique  which  converges  if  an  answer  is  possible  but  terminates  if  user-specified  limits 
are  reached.  The  run  time  is  dependent  upon  the  input  values  of  several  variables:  NC,  NREP, 
MMAX,  NMAX,  EPS,  and  EMIN.  The  user  will  learn  by  experience  how  to  choose  reasonable  input 
values.  The  example  problem  may  provide  some  insight. 

In  the  program  listing  of  figure  5,  the  routine  VRAND  calls  upon  the  routine  GGNOF, 
which  is  not  listed  because  it  is  a proprietary  program  from  the  International  Mathematical  and 
Statistical  Library.  The  properties  of  GGNOF  are  the  return  of  a random  normal  deviate  with  zero 
mean  and  unit  variance.  The  seed,  IS,  is  changed  with  each  call  of  the  routine.  The  function 
VRAND  converts  the  deviate  to  one  with  standard  deviation  SIG. 

The  input  value  of  SIG  should  be  the  root-mean-square  value  of  the  range  and 
deflection  standard  deviations;  i.e., 


SIG  range  + a~  deflection^  *4** 

This  choice  removes  the  effect  of  wind  direction  relative  to  the  line  of  flight. 

V.  SUMMARY. 

Methodologies  have  been  developed  and,  where  possible,  implemented  in  computer 
programs  for  estimating  smoke  munition  expenditure  requirements  under  a broad  range  of 
operational  conditions.  The  techniques  are  based  on  idealized  geometries  and  simplifying 
assumptions  which  are  expected  to  over-estimate  requirements  by  an  undeterminable  margin. 
Validation  of  the  estimates  is  desirable  but  has  not  been  attempted.  The  techniques  are  designed  to 
be  used  with  any  smoke-effect  simulation  model  capable  of  providing  information  on  effective 
screen  size  as  a function  of  time. 
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e : XQT  F.MUNX 
RANDOM  NUMBER  SEED 
INITIAL  PUINTS 
REPLICATIONS 
MAXIMUM  ROUNDS 
MAXIMUM  POINTS 
MAXIMUM  RELATIVE  ERROR 
SUBMUNIT1CNS  PER  ROUND 
RADIUS  OF  IMPACT  ZONE,  M. 

HALF-WIDTH  OF  SUBMUNITION  SCREEN,  M. 
LENGTH  OF  FRONT,  M. 

DELIVERY  ERROR(SIGMA) , M. 

MINIMUM  COVERAGE 


Figure  4.  MUNX  Program  Output 


1234567 

10 

30 

30 

20 

.500-03 

4 

74.00 
1 . 9C0 
125.00 
5.00 
.950 


INDS 

COVERAGE 

VARIANCE 

POINTS 

REL.  ERR. 

REP.  VAR. 

1 

.107+00 

.301-03 

16 

.243-03 

.965-06 

2 

.204+00 

.859-03 

16 

.144-03 

.525-06 

3 

.288+00 

.156-02 

16 

.141-03 

.186-05 

4 

.365+00 

.223-02 

16 

.143-03 

.741-06 

5 

.432+00 

.276-02 

16 

.141-03 

.148-05 

6 

.492+00 

•326-02 

16 

.143-03 

.154-05 

7 

.546+00 

.358-02 

16 

.135-03 

.129-05 

8 

.594+00 

. 382-02 

16 

.139-03 

.425-05 

9 

.637+00 

. 382-02 

16 

.127-03 

.180-05 

10 

.675+00 

.385-02 

16 

.128-03 

.191-05 

11 

.709+00 

.378-02 

16 

.125-03 

. 112-05 

12 

.739+00 

.370-02 

16 

.124-03 

.139-05 

13 

.766+00 

.349-02 

16 

.119-03 

.452-06 

14 

.791+00 

.328-02 

16 

.117-03 

.682-06 

15 

.812+00 

.307-02 

16 

.112-03 

.617-06 

16 

.832+00 

.285-02 

16 

.110-03 

.617-06 

17 

.849+00 

.259-02 

16 

.103-03 

.518*06 

18 

.864+00 

.241-02 

16 

. 102-03 

.551-06 

19 

.878+00 

.214-02 

16 

.983-04 

.286-06 

20 

.890+00 

.196-02 

16 

.942-04 

.469-06 

21 

.902+00 

.177-02 

16 

.932-04 

.370-06 

22 

.912+00 

.157-02 

16 

.894-04 

.345-06 

23 

.920+00 

.142-02 

16 

.846-04 

.286-06 

24 

.928+00 

.125-02 

16 

.809-04 

.271-06 

25 

.936+00 

. 112-02 

16 

.787-04 

.337-06 

26 

.942+00 

.983-03 

16 

.771-04 

.181-06 

27 

.948+00 

.863-03 

16 

.725-04 

.206-06 

28 

.953+00 

.754-03 

16 

.694-04 

.148-06 

32 
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Table  4.  MIJNX  Input 

Card  No. 

Variable 

Type 

Definition 

I 

1 

IS 

Integer 

Random  number  seed,  odd  integer. 

NC 

Integer 

Initial  number  of  points  for  midpoint 

integration  scheme. 

NREP 

Integer 

Number  of  replications. 

IfiAX 

Integer 

Maximum  number  of  rounds. 

NMAX 

Integer 

Maximum  value  of  NC. 

EPS 

Real 

Maximum  allowable  relative  error  in 

* 

integration. 

* 

2 

N 

Integer 

Number  of  submunitions  per  round. 

R 

Real 

Radius  of  impact  area  (one  round),  m. 

i 

W 

Real 

Half-width  of  one  submunition  screen,  m. 

L 

Real 

Length  of  front,  m. 

SIG 

Real 

Standard  deviation  of  delivery  system,  m. 

i 

EMIN 

Real 

Minimum  expected  fractional  coverage 

of  front. 
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11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 


.MUNX 

DIMENSION  V( 1000) 

REAL  L 

READ!  5, 100)  I S , NC  , NREP  , MivlAX , NMAX  , E PS 
1 RE AO (5, 100,END*99)  N , R , W , L , S IG , EMI N 

WRITE! 6, 101 ) IS.NC , NREP, MMAX.NMAX.EPS, N.R.W, L.SlG, EMIN 
WR I T E ( 6 , 1 03 ) 

L=.5*L 

DO  50  M=1 , MMAX 

EER=0 

SEE  = 0 

SE  = 0 

SV*0 

DO  40  NRP* 1 , NREP 
' DO  10  I x 1 , M 
10  V( I )=VRAND( IS.SIG) 

EEE  = 0 
NP=NC 
IT  EE=0 
VR  = 0 

DO  15  J=1 ,NP 
X=-l+(2*J-1 )*L/NP 
PR  = 1 

DO  20  K=  1 , M 
A*X-W-V(K) 

B=X+W-V(K) 

P=1-G( B.R)+G( A.R) 

PR  = P R*  P 

IFfPR.lE.O.)  GO  TO  21 
20  CONTINUE 
2 T EE  = EE  + PR**N 
15  VR=VR+PR**(2*N) 

EEF=1-EE/NP 

ERR=ABS( EEF-EEE)/( EEF+1 .E-20) 

IF(ERR.LT.EPS)  GO  TO  25 

EEE=EEF 

NP=NP+ 1 

IF(NP. LE.NMAX)  GO  TO  11 
25  NC=NP~1 
• EER=EER+ERR 
VR=VR/NP-( EEF-1 )**2 
SE=SE+EEF 
SEE=SEE+EEF**2 
40  SV=SV+VR 
E-SE/NREP 
VAR=  SV/NRE  P 
ER--EER/NREP 

RVAR  = (SEE-NREP*E**2)/(NREP-1  ) 

WRITEI6.102)  M.E.VAR.NP.ER.RVAR 
IF(E.GE.EMIN)  GO  TO  60 
50  CONTINUE 
60  GO  TO  1 
99  STOP 

100  FORMAT() 

102  FORMAT (5X. 15 , 2E1 0 . 3 , 5X , 1 5 , 2E 1 2 . 3 ) 

103  F0RMAT(//,5X' ROUNDS' ,T14, 'COVERAGE' ,T23, 'VARIANCE' ,T35 
* 'POINTS' ,T45, 'REL.  ERR . ' , T56 , ' REP . VAR.') 

101  FORMAT ( ' RANDOM  NUMBER  SEED ' , T60 , 1 1 0 , / , 


Figure  5.  MUNX  Program  Listing 


58 

* ' 

INITIAL  POINTS' , T60, 110,/. 

59 

* 1 

REPLICATIONS' ,T60, 110,/, 

60 

* ' 

MAXIMUM  ROUNDS' ,TG0, 110./, 

61 

* 1 

MAXIMUM  POINTS ', T60, I 1 0./, 

62 

♦ ' 

MAXIMUM  RELATIVE  ERROR ' , T60 . E 1 0 . 3 , / , 

63 

* ' 

SUBMUNITIONS  PER  ROUND '. T60 , I 1 0 ,/ , 

64 

* ' 

RADIUS  OF  IMPACT  ZONE,  M . ' , T60 , F 1 0 . 2 , / 

65 

* 1 

HALF-WIDTH  OF  SUBMUNITION  SCREEN,  M.‘, 

66 

* 1 

LENGTH  OF  TRONT,  M . ' , T60 , F 1 0 . 2 , / , 

67 

* 1 

DELIVERY  ERROR ( S I GMA ) , M . ' , T60  , F 1 0 . 2 , / 

68 

* 1 

MINIMUM  COVERAGE' .T60.F10. 3) 

69 

END 

SAD*  SMOKEX ( 1 ) .MUNXG 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


FUNCTION  G(X.R) 

IF(X.GE.R)  GO  TO  5 
IFiX.LE.-R)  GO  TO  10 
G-X‘SQRT(R**2-X**2 )+R»*2*ASIN(X/R) 
G=G/(3. 1415927*R**2) 

1  RETURN 
5 G=  .5 
GO  TO  1 
10  G*-.5 
GO  TO  1 
END 


SAD*  SMOKEX ( 1 ) . MUNX  R AND 

1 - FUNCTION  VRAND(IS.SIG) 

2 VRAND=SIG*GGNOF( IS) 

3 RETURN 

4 ENO 


Figure  5.  (Continued) 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS. 


Use  of  the  methodologies  and  computer  programs  in  this  report  has  not  produced  any 
estimates  which  appear  unreasonable  to  the  author,  who  has  witnessed  dynamic  test  firings  of  many 
smoke  munitions.  However,  there  is  no  familiarity  with  the  capabilities  of  fire  units  to  perform  in  a 
manner  compatible  with  the  methodology  assumptions.  For  this  reason,  the  only  recommendation 
to  be  made  is  that  analysts  with  the  required  expertise  should  evaluate  the  proposed  methodologies 
from  their  own  perspective  and  report  the  evaluation  in  the  literature. 
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